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INTRODUCTION: 


Only  a  relatively  small  amount  of  work  has  been  published  [1*3]  on 
the  transient  response  of  high  power  battery  systems  in  the  short  time 
domain  (O.l/is-lOms) .  Most  laboratory  studies  deal  with  scaled  down  versions 
of  large  cells  that  are  in  practice  usually  employed  to  provide  standby 
power  or  for  applications  in  which  continuous  operation  is  required,  rather 
than  for  pulse  generation.  In  tha  present  work,  the  impedance  of  an 
individual  elementary  cell  has  been  calculated  over  a  wide  range  of 
frequencies  in  order  to  show  the  effects  of  various  physical  and  geometric 
factors  that  are  significant  at  short  discharge  times. 


The  geometry  chosen  for  the  present  calculations  is  the  semi- 
infinite  parallel  plate  cell  with  infinitely  thick  solid  electrodes;  a 
configuration  commonly  referred  to  in  the  electronics  literature  as  a  strip- 
line.  Such  a  cell  is  illustrated  schematically  in  Fig.  1.  In  order  to  carry 
out  the  Impedance  calculation  this  battery  cell  has  been  modelled  as  a 
transmission  line  containing  uniformly  distributed  values  of  resistance, 
capacitance  and  inductance;  quantities  that  can  be  determined  directly  from 
the  physical  characteristics  chosen  for  the  model.  The  transmission  line 
can  be  analysed  as  a  distributed  network  of  differential  elements  [4] ,  each 
having  the  equivalent  circuit  shown  in  Fig.  2.  Its  electrical 
characteristics  can  then  be  expressed  in  terms  of  a  characteristic  impedance 
ZQ,  an  attenuation  constant  a,  and  a  phase  shift  (}. 
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where  r8  is  the  series  resistance,  gp  is  the  shunt  conductance,  Ls  is  the 
series  inductance,  and  Cp  is  the  shunt  capacitance,  all  per  unit  length. 


The  propagation  constant  7  for  the  transmission  line  shown  in  Fig. 
2  is  given  by  the  expression 


[  (rB  +  JwLg)  *  (gp  +  Jfc>Cp>  J1/2 
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Hw  propagation  constant  can  also  be'  written  in  terns  of  its  real  and 
imaginary  parts,  as 

T  -  a  +  j  fi  ....(3) 

The  real  part,  o,  is  the  attenuation  constant  and  Is  related  to  the 
penetration  length  of  the  perturbation  (along  the  z  axis  in  Figure  1).  fi  is 
the  phase  constant.  The  attentuation  constant  is  given  by  the  expression 

o  -  rs/2Zo  +  - (4) 

For  a  particular  frequency  all  but  1/e  of  the  total  current  is  drawn  from  a 
distance  1/a  fron  the  terminals.  This  distance  is  referred  to  as  the 
penetration  length.  The  ratio  of  currents  between  points  0  and  X,  separated 
by  a  distance  x,  for  a  wave  going  from  0  to  X,  is 

lx  /  Io  “  •*ox  ....(5) 

In  the  corresponding  time  domain,  this  implies  that  at  short  times  almost 
all  the  current  is  generated  within  about  one  penetration  length  (i.e.,  in 
the  regions  closest  to  the  current  collection  terminals) .  Regions  further 
along  the  z  axis  in  Fig.  1  make  little  or  no  contribution.  The  utilization 
Increases  with  tine ,  but  for  short  times  a  cell  need  have  a  length  of  only 
on  the  order  of  2  or  3  times  the  penetration  length  in  order  to  show  the 
behaviour  of  a  semi- ini fite  cell. 

The  above  treatment  yields  the  Impedance  as  a  function  of 
frequency,  but  what  is  preferred  in  practice  is  the  current/voltage 
response  as  a  function  of  tine  when  a  cell  or  battery  is  placed  on  load.  In 
principle  it  should  be  possible  to  compute  the  transient  response  by 
multiplying  the  impedance  function  by  the  transform  of  the  perturbation  and 
carrying  out  an  inverse  transformation  on  the  product.  This  will  be  the 
subject  of  a  future  paper. 


ESTIMATION  OF  PARAMETERS 

There  are  four  parameters,  rg,  gn,  L,  and  Cp  in  the  lumped  network 
representation  of  a  unit  section  of  transmission  line  shown  in  Fig.  1.  The 
physical  components  of  such  a  cell  are  Illustrated  schematically  in  Fig.  1, 
and  the  corresponding  equivalent  circuit  is  given  in  Fig.  3.  In  order  to 
calculate  the  impedance,  Fig.  3,  the  equivalent  circuit  representation  of 
the  physical  components,  must  be  translated  into  the  form  of  Fig.  1. 

The  values  of  the  parameters  in  Fig.  3  can  be  obtained  from  the 
physical  properties  chosen  for  the  model  as  follows:  rg  and  L,  are  modified 
by  the  skin  effect  [5].  The  magnetic  field  resulting  from  current  flow  in  a 
cell  configured  as  a  transmission  line  confines  the  current  to  the  outer 
layer  (skin)  of  the  conductors.  This  is  seen  at  high  frequencies  as  an 
exponential  decrease  in  current  density  in  a  conductor  as  the  distance  a way 
from  the  surface  increases.  The  equivalent  skin  depth,  d,  is  defined  as  that 
depth  over  which  a  current  density  equal  to  0.707  of  the  surface  current 
density  can  be  considered  to  be  distributed  uniformly. 
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4  -  V(<wr£#»)V2 

viwr*  a  la  tha  conductivity  In  Q'^.a'3-,  f  is  the  frequency  in  Hz,  and  n  is 
the  permeability  in  Henries. a" 1.  Thus,  even  for  an  infinitely  thick  plate, 
at  high  frequencies  all  but  1/e  of  the  current  will  be  carried  by  a  depth  d 
of  the  outer  region.* 

The  series  inductance  arises  froa  the  geometry  of  the  conductors 
and  is  given,  when  the  skin  depth  d  (or  the  thickness  of  the  electrodes 
theaselves  for  finite  thickness  electrodes)  is  saall  with  respect  to  the 
cell  width  V,  by 

L,  -  1.255xlO‘6*u/W  _ (7) 

where  u  is  the  electrode  spacing,  and  Lg  has  units  of  Henries/a.  The  series 
resistance  for  the  Infinitely  thick  electrodes  assuaed  is  dominated  by  the 
skin  effect  resistance  rsfc,  which  is  given  by 

r,k  -  3.974xlO*3-pmV2.fl/2/w  _ (8) 

where  is  the  resistivity  of  both  electrodes  in  Q-a,  and  rsk  has  units  of 
0/a.  In  the  present  model  it  is  assumed  that  all  plates  and  current 
collectors  are  Infinitely  thick,  so  rs  is  determined  solely  by  rs^  which  in 
turn  depends  on  the  skin  depth  and  on  the  conductivity  of  the  electrodes 
(i.e.  the  calculations  assunw  the  best  case).  For  example,  the  skin  depth 
for  copper  at  6xl05Hz  is  0.6mm  and  the  resistance  increases  above  that 
calculated  using  equation  (7)  as  the  plate  thickness  approaches  and  becomes 
saaller  than  this  value.  It  should  be  noted  that  for  electrodes  of  lower 
conductivity  (because  of  either  intrinsic  CQnductivity  or  porosity)  the 
skin  depth  increases  in  proportion  to  the  square  root  of  the  resistivity.  In 
consequence,  the  effective  resistance  will  only  increase  in  proportion  to 
the  square  root  of  an  increase  in  electrode  resistivity. 

It  should  be  noted  that  rsjc  is  a  function  of  frequency  although 
by  definition  a  real  resistance  is  frequency  invariant,  yet  the  phase  angle 
reaalns  zero.  This  is  in  apparent  contradiction  of  the  Kraaers-Kroenlg 
relations  [7],  but  it  must  be  recognised  that  rsk  is  only  an  effective 
resistance  and  is  a  result  of  the  nonunifora  distribution  of  current  density 
into  a  conductor  at  AC  frequencies. 

The  parallel  components  comprise  terms  for  the  electrode 
Interfaces  in  series  with  terms  for  the  electrolyte.  The  interfacial 
components  are  represented  by  the  parallel  combination  of  a  Faradaic 
resistance,  rp,  and  a  double-layer  capacitance,  C<jx.  The  electrolyte  terms 
are  the  Ohmic  resistance,  rg,  and  in  parallel  the  dielectric  capacitance, 
Cg,  where  [5] 

Cg  -  «.«e. W/u  - (9) 


♦NOTE :  The  skin  depth  is  often  surprisingly  small.  Sixty  hertz  power  I 
transmission  lines  are  rarely  larger  than  two  skin  depths  in  diaaeter  1 
(d£u(60Hz]  “  6am).  Increasing  the  diaaeter  of  the  wire  decreases  the  I 
resistance  in  proportion  to  the  circumference  rather  than  to  the  cross-  | 
sectional  area. [6]  3 
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«  is  the  dielectric  constant  of  the  electrolyte,  and  C£  is  in  units  of 
Farads/a?.  In  the  nodal  discussed  here,  the  interface  is  considered  to  be 
saooth. 

The  series  inpedances  for  the  interfaces  and  for  the  electrolyte 
nust  then  be  combined  in  series  and  then  converted  to  a  coaplex  admittance 
of  the  fora  gp  ♦  JwCp  for  use  in  equation  (3).  It  should  be  noted  that  gp 
and  Cp  here  do  not  correspond  to  any  siaple  physical  quantities. 


RESULTS  AND  DISCUSSION  : 

The  strip-line  cell  shorn  in  Fig.  1  has  been  modelled  as  a 
transmission  line.  The  impedance  per  unit  width  and  length  has  been  computed 
as  a  function  of  frequency  for  different  values  of  electrode  and  electrolyte 
conductivity  and  cell  thickness.  For  the  purposes  of  these  initial 
calculations,  the  case  where  no  Faradaic  processes  are  occurring  has  been 
considered;  that  is  the  Faradaic  resistance,  rp  in  Fig.  3,  has  been 
considered  to  be  infinite.  A  constant,  frequency- independent  double- layer 
capacitance  of  50uF/cm^  has  been  used,  and  it  has  been  assumed  that  the 
electrodes  are  infinitely  thick  so  that  the  series  resistance  is  equal  to 
r«k-  The  cell  is  assumed  to  be  Infinitely  long,  although  only  a  small 
fraction  of  this  length  will  actually  be  able  to  deliver  current  to  a  load 
at  short  times  or  high  frequencies.  The  penetration  length  has  been 
calculated  as  a  function  of  frequency  for  each  set  of  conditions. 

The  results  are  presented  as  Bode  plots  for  log  Z,  the  logarithm 
of  the  modulus  of  the  coaplex  impedance  (equivalent  to  log  |Z0|),  the  phase 
angle  i ,  and  the  penetration  length  (equal  to  1/a  -  equation  (5)). 

Depending  on  the  values  of  the  physical  variables  chosen  for 
the  calculation,  the  equivalent  circuit  shown  in  Figure  2  can  be  simplified 
to  yield  a  number  of  limiting  cases.  If  wLg  and  wCp  are  small  with  respect 
to  rs  end  g«,  respectively,  then  the  characteristic  Impedance,  Z0,  will  be 
resistive  in  nature.  If  wLs  and  uCp  are  large  with  respect  to  rs  and  gp, 
then  the  circuit  will  reduce  to  an  LC  network  and  Z0  will  again  be 
resistive.  If  uC  is  small  with  respect  to  gp  only  then  the  circuit  reduces 
to  an  RL  network  and  ZQ  will  show  a  positive  phase  shift,  while  if  wL  is 
small  with  respect  to  rs  only  then  the  circuit  reduces  to  an  RC  network  and 
Z0  has  a  negative  phase  shift.  The  results  of  the  calculations  presented 
here  indicate  that,  depending  on  the  frequency  range  of  interest  and  the 
values  of  the  physical  variables  chosen,  a  number  of  these  different 
limiting  cases  as  well  as  intermediate  behaviour  will  be  observed  with  real 
cells. 

Figure  4  is  a  set  of  Bode  plots  calculated  for  electrode 
resistivities  ranging  from  1.7x10"®  to  1.7X10"1  O.cm.  It  may  be  noted  that 
Cu  has  a  resistivity  of  1.7xl0‘®  O.cm.  The  cell  thickness  is  0.01cm,  while 
the  electrolyte  resistivity  has  been  set  high  at  l.OxlolO  O.cm.  This  latter 
value  gives  the  limiting  case  where  wCp  is  large  with  respect  to  gp.  Figure 
4  shows  that  at  high  frequencies,  for  the  lower  values  of  pf wL8  is  much 
greater  than  the  series  resistance,  r8,  and  the  behaviour  of  the  cell 
approaches  that  of  a  pure  LC  strip  line.  (It  can  be  seen  from  equation  (1) 
that  if  rB-gp-0 ,  then  the  jw  factors  cancel  and  ZQ  becomes  wholly  real.)  Z 


approaches  a  constant,  purely  resistive,  value  of  about  50  (determined 
largely  by  the  values  chosen  for  the  cell  thickness  and  the  electrolyte 
dielectric  constant),  while  9  tends  to  zero. 

For  more  resistive  electrodes  and  at  lower  frequencies,  the 
contribution  from  the  series  resistance  term  is  larger,  and  the  overall 
behaviour  becomes  capacltative .  The  impedance  increases,  while  9  tends  to  • 
45°.  For  Cu,  for  example,  Z  begins  to  increase  at  a  frequency  of  about 
lO^Hz.  It  might  be  expected  that  as  rs  becomes  dominant  at  low  frequencies 
(since  rgfc  decreases  in  proportion  to  only  the  square  root  of  f)  and/or  as 
the  resistivity  of  the  electrode  increases,  the  behaviour  of  the 
transmission  line  would  approach  that  of  an  RC  network.  For  this  case  the 
Kramers -Kroenig  relations  [4]  predict  that  when  9  goes  to  -45°,  log  Z  should 
decrease  with  frequency  with  a  slope  of  -1/2.  However  Figure  4  shows  that 
the  computed  value  of  log  Z  only  decreases  with  a  slope  of  -1/4.  The 
probable  explanation  for  this  behaviour  is  that  not  only  the  values  of  the 
circuit  components  but  also  the  equivalent  circuit  representation  that  is 
applicable  changes  with  frequency.  This  result  provides  a  further 
illustration  of  the  necessity  of  predicting  the  time  or  frequency  dependent 
behaviour  of  distributed  systems  from  computations  based  on  a  complete 
model,  rather  than  on  intuitive  extrapolations  from  steady  state  behaviour. 
It  is  not  in  general  valid  to  replace  a  distributed  network  by  a  single 
simple  equivalent  circuit.  In  the  model  used  here,  r8  continues  to  change 
with  frequency  because  an  infinite  electrode  thickness  has  been  assumed.  If 
instead,  a  finite  plate  thickness  is  assumed,  then  once  the  skin  depth 
exceeds  this  thickness  (at  low  frequencies  -  see  equation  (8)),  rs  should 
tend  to  a  constant  value  and  the  slope  of  log  Z  should  tend  to  -1/2. 

The  penetration  lengths  are  also  shown  in  Fig.  4,  and  it  can  be 
seen  that  for  this  case  of  a  highly  resistive  electrolyte  they  are 
relatively  large  over  almost  the  whole  frequency  range.  At  low  frequencies 
they  decrease  by  less  than  an  order  of  magnitude  as  the  electrode 
resistivity  increases  by  six  orders  of  magnitude.  At  the  higher  frequencies 
electrode  resistivity  has  almost  no  effect  on  the  penetration  length.* 

Figure  5  shows  the  same  calculation  as  that  shown  in  Figure  4 
performed  for  a  cell  thickness  of  0.1  cm.  It  can  be  seen  that  the  resulting 
set  of  Bode  plots  have  the  same  shape ,  but  that  the  impedance  curves  have 
been  shifted  upwards  (to  higher  impedances),  and  both  the  log  Z  and  the  9 
curves  have  been  translated  to  the  left,  i.e.  to  lower  frequencies.  The 
penetration  length  is  increased  slightly  at  lower  frequencies.  These  changes 
in  the  limiting  behaviour  are  the  result  of  the  order  of  magnitude  decrease 
in  the  dielectric  capacitance  of  the  cell  corresponding  to  the  order  of 
magnitude  increase  in  electrode  spacing.  Because  this  quantity  is  much 
smaller  than  the  double-layer  capacity,  it  dominates  the  value  of  Cp.  (See 
Figure  3.)  The  limiting  value  of  Z  increases  because  Cp  decreases  and  also 
because  Ls  increases  with  the  increase  in  spacing  (equation  (1)). 


*  NOTE:  The  calculation  above  corresponds  to  the  case  of  a  cell  having  a 
hlghly-resistive  (e.g.  nonaqueous)  electrolyte  and  indicates  that 
significant  currents  can  still  be  drawn.  The  calculation  also  applies  to  a 
capactitor  having  a  configuration  corresponding  to  this  strip  line,  which 
will  also  demonstrate  resistive  behaviour  at  high  frequencies  and  a 
significant  increase  in  impedance  at  low  frequencies. 


Figure  6  shows  the  effect  of  varying  the  electrolyte  resistivity 
over  the  range  10*2  to  103  n.cm  when  the  electrodes  are  highly  conductive. 
The  cell  thickness  is  0.1  cm,  and  p^  is  1.7x10*®  n.cm.  Over  a  wide  frequency 
range,  from  10^  to  10®  Hz  depending  on  the  value  of  pg,  gp  is  much  larger 
than  wCp  and  the  transmission  line  shows  a  positive  phase  shift.  This  region 
is  characterised  by  log  Z  increasing  with  increasing  frequency  with  a  slope 
of  1/2,  while  the  phase  angle  tends  to  +45°.  This  is  the  behaviour  expected 
for  an  LR  stripline.  The  frequency  at  which  this  behavior  is  seen  depends  on 
the  value  of  pg,  shifting  to  lower  frequencies  as  pg  increases.  At  still 
lower  frequencies  r8  becomes  large  with  respect  to  wLs,  and  log  Z  reaches  a 
constant  value  (about  0.002  Q)  characteristic  of  purely  resistive  behaviour. 
Correspondingly,  9  decreases  to  zero. 


At  frequencies  above  the  LR  region  wCp  becomes  increasingly 
significant,  and  as  it  exceeds  gp  the  behaviour  of  the  line  again  becomes 
resistive  (see  the  discussion  of  an  LC  network  above).  Log  Z  tends  to  a  new 
limiting  value,  almost  two  orders  of  magnitude  greater,  and  9  decreases  to 
zero. 


Figure  6  also  shows  that  the  penetration  length  is  strongly 
dependent  on  the  electrolyte  resistivity,  decreasing  as  pg  decreases. 


Figure  7  shows  the  effect  of  electrolyte  resistivity  when  the 
electrodes  have  a  relatively  low  conductivity.  Again  pg  ranges  from  10*2  to 
10®  0. cm,  while  the  electrode  resistivity  is  3.5xl0*3  n-cm,  corresponding  to 
an  electrode  material  such  as  carbon.  The  electrolyte  thickness  is  0.1  cm. 

Over  limited  frequency  ranges,  the  behavior  of  the  cell  can  be 
described  by  simplified  versions  of  the  equivalent  circuit  shown  in  Fig.  2, 
and  these  are  illustrated  in  Fig.  7  for  the  curve  corresponding  to  an 
electrolyte  resistivity  of  10^  tl-cm.  At  high  frequencies  (region  E) ,  this 
transmission  line  behaves  like  an  LC  network  and  ZQ  has  the  characteristics 
of  a  pure  resistance  (see  equation  (1)).  Log  Z  is  constant,  and  the  phase 
shift,  9,  approaches  zero.  As  the  frequency  decreases  (region  D)  wCp  becomes 
less  than  gp  and  the  behavior  becomes  that  of  the  LR  circuit  shown.  The 
impedance  decreases,  and  9  becomes  positive.  As  the  frequency  continues  to 
decrease  (region  C)  wLs  becomes  comparable  to  rs  and  the  equivalent  circuit 
tends  to  become  purely  resistive  as  shown  (region  B) .  At  the  lowest 
frequencies  shown  (region  A)  the  interfacial  terms  begin  to  become 
important.  Here  starts  to  become  the  dominant  parallel  component  and 
the  Impedance  increases  again,  while  9  becomes  negative.  This  may  be  a 
consequence  of  the  assumption  of  infinitely  long  plates  and  infinitely  thick 
electrodes.  As  the  skin  depth  continues  to  increase  with  decreasing 
frequency,  the  penetration  length  may  become  limiting.  It  should  be  noted 
that  the  Kramers -Kroenig  relationship  does  not  hold  here  either,  since  a 
depends  on  r8  which  is  frequency  dependent. 

As  in  the  previous  example  the  penetration  length  is  dominated, 
particularly  at  high  frequencies,  by  pg.  The  effect  of  the  lower  value  for 
Plf  appears  as  a  slight  decrease  in  the  penetration  length  at  low 
frequencies . 


From 


the  examples  above,  it  may  be  concluded  that  at  high 


frequencies /short  tines  the  behaviour  of  the  simple  cell  modelled  here  is 
United  by  the  impedance  resulting  froa  the  physical  geometry  of  the  cell 
components.  These  factors  must  therefore  be  given  strong  consideration  when 
designing  batteries  for  high-power  short -duration  pulse  applications.  The 
continuation  of  this  work  is  now  in  progress  to  extend  the  model  to  take 
into  account  the  effects  of  the  electrochemical  interface  and  more  complex 
cell  geometries. 


SUMMARY  : 

1.  The  impedance  resulting  from  the  physical  geometry  of  the  components 
will  limit  the  high  frequency/short  time  performance  of  a  cell  and  must 
therefore  be  considered  in  addition  to  the  electrode  kinetics  in  designing 
short  duration  pulse  batteries. 

2.  In  order  to  take  these  effects  into  account  in  estimating  battery 
performance,  it  is  necessary  to  model  the  battery  as  a  complex  distributed 
network.  The  behavior  of  these  networks  does  not  follow  intuitive  concepts 
applicable  to  steady- state  behaviour  or  to  systems  that  can  be  uniquely 
represented  by  a  simple  equivalent  circuit.  Depending  on  the  condition 
chosen,  any  of  the  variables  considered,  i.e.,  electrode  resistance, 
electrolyte  resistance,  or  cell  thickness,  can  affect  the  performance  of  a 
cell. 

3.  The  calculated  values  of  penetration  length  indicate  that  at  short 
times/high  frequencies  current  can  only  be  drawn  from  the  regions  closest  to 
the  terminals  of  most  conventional  large  battery  cells,  so  that  only  a 
fraction  of  their  total  power  is  available.  This  problem  becomes  greater  as 
the  size  of  the  battery  is  increased  since  the  ratio  of  1/a  to  cell 
dimensions  becomes  increasingly  unfavorable.  Transients  measured  for  small 
test  cells  should  not  necessarily  be  expected  to  scale  as  the  cell/battery 
size  is  increased. 
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FIGURE  2  :  Equivalent  Circuit  of  Individual  Element  in  Distributed 

~~  Network  Analysis  of  Strip-line  Cell 


FIGURE  3  Equivalent  Circuit  Representation  of  Physical  Components 

of  Strip-line  Battery  Cell 
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